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ABSTRACT 

The enormous theoretical potential of Quantum Information Processing (QIP) is driving the 
pursuit for its practical realization by various physical techniques. Currently Nuclear Magnetic 
Resonance (NMR) has been the forerunner by demonstrating a majority of quantum algorithms. In 
NMR, spin systems consisting of coupled nuclear spins are utilized as qubits. In order to carry out 
QIP, a spin system has to meet two major requirements: (i) qubit addressability and (ii) mutual 
coupling among the qubits. It has been demonstrated that the magnitude of the mutual coupling 
among qubits can be increased by orienting the spin-systems in a liquid crystal matrix and utilizing 
the residual dipolar couplings. While utilizing residual dipolar couplings may be useful to increase 
the number of qubits, nuclei of same species (homonuclei) might become strongly coupled. In 
strongly coupled spin-systems, spins loose their individual identity of being qubits. We propose 
that even such strongly coupled spin-systems can be used for QIP and the qubit-manipulation can 
be achieved by transition-selective pulses. We demonstrate experimental preparation of pseudopure 
states, creation of maximally entangled states, implementation logic gates and implementation of 
Deutsch-Jozsa (DJ) algorithm in strongly coupled 2,3 and 4 spin systems. The energy levels of the 
strongly coupled 3 and 4 spin systems were obtained by using a Z-COSY experiment. 



I. INTRODUCTION 



The theoretical success of exploiting the quantum nature of physical systems in certain information 
processing tasks like prime factorization [1] and unsorted database search [2] has motivated the 
pursuit for the practical realization of Quantum Information Processing (QIP) [3-5]. With the 
demonstration of many quantum algorithms, Nuclear Magnetic Resonance (NMR) is now considered 
as a suitable test-bed for QIP. One of the main challenges for the progress of NMR QIP is "how 
to increase the number of qubits?". In this direction several attempts are being made, such as (i) 
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find molecules with different chemical shifts and J-couplings, and (ii) use of dipolar and quadrupolar 
couplings. This paper concentrates on one aspect: "how to use dipolar couplings among homonuclear 
spins?" . This problem is outlined in the following paragraphs. 

In NMR, systems consisting of coupled spin- 1/2 nuclei form qubits. In order to carry out QIP, the 
spin-system has to meet two main requirements: qubit addressability and mutual coupling among 
the qubits. In liquid state NMR using isotropic fluids, the qubit addressability is normally provided 
by the differences in Larmor frequencies of the various spin-1/2 nuclei, while the mutual coupling is 
normally provided by the scalar (J) coupling among the nuclei connected by covalent bonds. The 
Hamiltonian for a J-coupled spin system is [6], 

H = Hz + Hj 

^ ' bJjljz ^ ' 2,77 Jij {^lizljz lixljx ~\~ liyliy) (1) 

where Hz is the Zeeman Hamiltonian, and Hj is the coupling Hamiltonian. When 2nJij <C \uJi~ u>j\, 
the system is said to be weakly coupled, and the Hamiltonian can be approximated to [6], 

n = ^UiI iz + £ 27rJ ij /«/ ja (2) 

i i,j(i<j) 

For qubit addressability, all uii should be sufficiently dispersed and all Jij should be non-negligible 
(> 1Hz) and unequal in magnitude. In such a circumstance each spin can be treated as a qubit 
and the coupled nuclei as several qubits. The values of depend on the covalent bonds connecting 
spins i and j, and normally has a small range (< 10 2 Hz) and becomes too small (< 1 Hz) if the spins 
are connected by more than 4-5 covalent bonds. This places a natural limit on the number of qubits 
reachable by liquid state NMR using J-couplings alone. To overcome this limitation the possibility 
of using dipolar couplings was considered. The truncated Hamiltonian for dipolar interaction is [6,7], 

Kd= 2ttA,(3/^ - Ii.Ij) (3) 

i,j(i<j) 

The dipolar coupling D^, between spins of gyromagnetic ratios ji and 7^ whose inter-distance vector 
Tij makes an angle 8{j with the Zeeman magnetic field is of the form [6,7], 

A, = |^(l-3co S 2 %) (4) 

Dipolar couplings among common nuclear species have larger magnitudes (~ 10 3 Hz) and longer 
range than the scalar couplings. However, in isotropic liquids the time average of vanishes, while 
in solids there are too many dipolar couplings resulting in broad unresolved lines and loss of qubit 



addressability. In molecules oriented in a liquid crystal matrix, while the intermolecular dipolar 
couplings are vanishingly small, the intra molecular dipolar couplings survive, scaled down by the 
order parameter (SV,) of the liquid crystal [7,8], 

U iJ - 47rr 3 <! dCO8 v >- (5) 

In such systems one obtains a finite number of sharp well resolved spectral lines making it possible 
to use such systems for NMR-QIP. In the NMR-QIP experiments implemented so far, the systems 
have been chosen such that either (i) 27r(Jj J - + 2Dij) <C \(u>i — yielding weakly coupled spin- 
systems which is the case for the heteronuclear spin-systems [9,10] or (ii) the coupling 2(7^ + 2/}^) is 
finite and \u)i — ujj\ = 0, i.e., equivalent-spins case [11,12]. In the latter case, the symmetry filtering 
of energy levels become increasingly difficult for higher number of qubits [12]. Even though the 
heteronuclear spins oriented in liquid crystal matrix are excellent for QIP since they provide very 
good qubit addressability as well as large mutual coupling, the use of more than 3 to 4 heteronuclear 
spins is limited by the extensive hardware requirements. Therefore for reaching larger number of 
qubits, one needs to utilize homonuclear (nuclei of same species having same 7 but different chemical 
shifts ) spins oriented in a liquid crystal. 

Homonuclear spins oriented in a liquid crystal generally become strongly coupled since the dipolar 
couplings become comparable to or more than the differences in Larmor frequencies \u>i — In 
such a situation, the Zeeman and the coupling parts of the Hamiltonian do not commute. Therefore 
the eigenstates of strongly coupled spins are obtained as the linear combinations of product states of 
various spins and the individual spins can no more be treated as qubits. We propose and demonstrate 
here that, the 2 N eigenstates of a coupled N-spin 1/2 system can be treated as an N qubit system 
even in the presence of strong coupling. Similar idea has already been used in demonstration of QIP 
using quadrupolar (S >l/2) nuclei oriented in high magnetic field, where the 2S+1 non-equidistant 
energy levels have been treated as N-qubit systems, where 2 N = 2S + 1. So far S=3/2 and 7/2 have 
been utilized respectively as 2 and 3 qubit systems [14-18], for various NMR-QIP. 

While substantial work has been carried out in NMR-QIP using weekly coupled spin systems 
[19], till now the use of the strongly coupled spin systems for QIP has not been experimentally 
demonstrated, presumably because (i) spin-selective pulses are not defined in the case of strongly 
coupled spins [6] , and (ii) the difficulty in constructing a general unitary operator using the evolution 
under scalar coupling [20]. The problem of using scalar coupling evolution of a strongly coupled two- 
spin system for a general unitary transform has recently been addressed theoretically, but extending 



to a N-spin system is complicated [21]. However, we note that unlike spin-selective pulses, the 
transition-selective pulses are well defined even in strongly coupled spin-systems [6] and hence it 
is possible to construct a unitary transformation using transition selective pulses. In section II we 
demonstrate NMR-QIP on a strongly coupled two spin-system in isotropic medium by preparing 
pseudopure states, implementing DJ algorithm, by creating Einstein-Podolsky-Rosen (EPR) state 
and by implementing logic gates. Section III describes the creation of Greenberger-Horne-Zeilinger 
(GHZ) states and implementation of two-qubit DJ algorithm on a strongly coupled three spin system 
in an oriented medium after labeling the transitions using the Z-COSY experiment. Labeling of 
transitions, preparation of pseudopure states, and implementation of gates on a four-spin strongly 
coupled system are demonstrated in section IV. 
In this study the strongly coupled systems used are: 



SI. No. 


Sample 


Solvent 


Isotropic/ Oriented 


No. of qubits 


1 


Trisodium citrate 


D 2 


Isotropic 


2 


2 


Organometallic compound (I) 


CDC1 3 


Isotropic 


2 


3 


l-bromo-2,3-dichlorobenzene 


ZLI-1132 


Oriented 


3 


4 


2-chloroiodobenzene 


ZLI-1132 


Oriented 


4 



The experiments have been carried out on a Bruker DRX-500 NMR spectrometer at 300K temper- 
ature. 

II. TWO-SPIN STRONGLY COUPLED SYSTEM 

The four eigenstates of a strongly coupled two-spin system (of spin 1/2 nuclei; AB spin-system) in 
isotropic medium are, \aa), cos ©la/3) + sin 0|/3a), cosQ\/3a) — sin 0| a/5), and \f3f3) (Figure 1(a)), 
where O = \ iwr 1 3 ab I a — ^b)) [6]. These eigenstates are labeled respectively as |00), 1 01) , 
1 10) , and 1 1 1 ) , thus forming a two-qubit system (Figure 1(b)). To demonstrate QIP on such a system, 
we have taken the strongly coupled X H spins of trisodium citrate (Figure 1(c)). In this system, the 
scalar coupling (J) is 15 Hz, the difference in Larmor frequencies (A/) is 55.5 Hz, and the strong 
coupling parameter (O) is 7.6°. The equilibrium spectrum of the system is shown Figure 1(d). 

A. Preparation of Pseudopure states 

In QIP, the computation normally begins from a definite initial state known as a pure state [3,4]. 
In NMR however because of the small energy gaps, it is not possible to realize a pure state, wherein 



the whole population is in one energy level, since it requires very low temperatures as well as very 
high magnetic fields. However, an alternate solution was discovered to overcome this problem [22,23]. 
In thermal equilibrium, NMR density matrix can be written as 

p = 2~ N {I + ep dev }. (6) 

The first part is a normalized unit matrix which corresponds to a uniform population background. 
The second part containing the traceless deviation density matrix pd ev (with a small coefficient 
e ~ 10~ 5 ) evolves under various NMR Hamiltonians, and gives measurable signal. It was observed 
independently by Cory et al [22] and Chuang et al [23], that by applying certain pulse sequence to 
the system in equilibrium, we can prepare the so called pseudopure density matrix, 

Ppps = 2~ N {(1 - e'/2 N )I + e'p pure }. (7) 

The first part is again a scaled unit matrix, but the second part corresponds to a pure state. Such 
pseudopure states mimic pure states [22,23]. Many methods have been proposed for the preparation 
of pseudopure states including spatial averaging [22,24], temporal averaging [25], logical labeling 
[23,26,27], and spatially averaged logical labeling [28]. Some of the other methods include the 
preparation of pseudopure states via cat states [29] and preparation of pair of pseudopure states 
[11]. Both J-evolution and transition selective pulse methods have been utilized for preparation of 
pseudopure states [22,25,13,27,16,30]. 

We have adopted the method of spatial averaging using transition selective pulses. The Boltzmann 
distribution of populations at high-temperature approximation is linear with energy gap. The equi- 
librium populations (in excess of a large uniform background population) of a homonuclear two spin 
system is given in Fig 2(a). For creating a pseudopure state all the populations except one of the 
states have to be equalized. This distribution can be achieved by a sequence of transition selective 
pulses intermittent with field gradient pulses to destroy any coherence created in the process. 

A transition-selective pulse of nutation angle 9 and of any transverse phase between states (i, j), 
changes the populations as follows: 

K=^cos 2 (^/2)+^sin 2 (^/2) 

p'j = p 3 cos 2 (9/2) + Pi sin 2 (6/2). (8) 

To prepare the 1 00) pseudopure state (Fig 2(b)) from equilibrium (Fig 2(a)), we use a sequence 
[(0)|io>Hii) _Gf z -(90 o )l° 1 > < - > l 11 > -G' z ] (pulses are applied from left to right), and it is inferred from the 
deviation populations of Fig 2(a) and (b), that 9 should be such that cos 2 (6'/2)=2/3. This yields 9 = 



70.5°. 6 = 90° on the other hand equalizes the populations of the two levels to an average value. Thus, 
the |00) pseudopure state is prepared by the pulse sequence [(70.5°) |10> ~ |n> - G z - (90°) |01 > m|11 > - G' z ] 
and the corresponding spectrum is given in Fig. 2(b). The 1 01) and 1 10) pseudopure states (Figure 
2(c), 2(d)) are respectively prepared by applying respectively (180°) (|00) ^ |01)) and (180°) (|00) ^ |10)) 
pulse after creating the 1 00) pseudopure state. The 1 1 1) pseudopure state (Figure 2(e)) is prepared 
by the pulse sequence (70.5°)(l 00 >^l 01 > - G z - (90°) l°°>^l 10 > - G' z . The observed intensities in the 
spectra on the left hand side of Fig. 2 correspond to the created population distribution and hence 
confirm the creation of the pseudopure states. 

B. Deutsch-Jozsa Algorithm 

Deutsch-Jozsa (DJ) algorithm is one of the first quantum algorithms which successfully demon- 
strated the power of QIP [32,33]. The task of DJ algorithm is to distinguish between two classes 
of many-input-one-output functions, constant and balanced. Constant functions are those functions 
in which all the outputs are same independent of inputs; and balanced functions are those in which 
half the number of inputs give one output and the other half gives another output. Classically, given 
a function of n input bits, it takes 2 n ~ 1 + 1 function-calls on an average, to determine whether the 
function is constant or balanced, whereas DJ algorithm needs only one function-call for any number 
of qubits. DJ algorithm has been implemented in NMR using scalar coupling evolution as well as 
using spin and transition selective pulses [27,34-36]. We have followed Cleve's version of DJ algo- 
rithm which requires one extra work qubit [33]. The circuit diagram and the NMR pulse sequence for 
implementing 1-qubit DJ algorithm are shown in Figures 3(a) and 3(b) respectively. The experiment 
begins with the 1 00) pseudopure state. An initial (7r/2)_ H pulse (the pseudopure Hadamard operation 
[35]) on all the qubits creates a superposition. It may be noted that unlike weakly coupled spins the 
superposition created here is not uniform in the eigenbasis, since the coefficients of various eigen- 
states are different. However, as is shown here it is still possible to distinguish between the different 
classes of functions. The Hadamard operation is followed by an unitary operator U / corresponding 
to the given function /. The unitary operator carries out the transformation |r, s) — * \f(s) © r, s), 
where \r) and \s) are the states of the work-qubit and the input-qubit respectively. The four different 
one-qubit functions and corresponding unitary operator as well as r.f. pulses are listed in Table 1. 
The test for a balanced function is that the transitions of input-qubit will gain opposite phases at 
the end of the algorithm. The experimental results corresponding to all the four functions /i, /2, 



/ 3 , and along with their corresponding simulated spectra are given in Figures 3(c)-3(f). From the 
spectra we can identify that functions f\ and f 2 are constant since the transitions of 1 and 2 are of 
same phase, whereas / 3 and / 4 are balanced since the transitions of 1 and 2 are of opposite phase. 



C. Creation of an EPR state 



Einstein-Podolsky-Rosen (EPR) pairs are the maximally entangled pairs of the form 

-L(|00>±|11» or 



V2 



(|01)±|10» 



(9) 



which are not reducible into product states of individual qubits [3]. The non-local correlation ex- 
hibited by these pairs have no classical equivalence, and are exploited in many branches of quantum 
information processing including quantum computation and quantum teleportation [3,4]. EPR states 
of a pair of weakly coupled nuclear spins have been earlier created by NMR using spin selective pulses 
and evolution of coupling [9]. 

Here we demonstrate the creation of EPR state on the above strongly coupled two-spin system 
using transition selective pulses and tomograph the result using non-selective pulses. Starting from 
1 00) pseudopure state, the EPR state ( 1 00) + \ ll))/y/2 can be created by applying the pulse sequence 
(pulses are to be applied from left to right) 



7T\ k 



(10) 



where k and I are the transitions, and <pi and <p 2 are the phases as shown in Table 2. For example, 



the unitary operator U for the pulse sequence 



is 



U = exp (iTT/WHii)) ■ exp (-i|4°°>~' 10 >) 



/ 



1 



1 -% 

V2 

iV2 

1 i 



where the operators are however applied from right to left. EPR state is obtained by applying U on 
1 00) pseudopure state 
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A phase-cycle over different combinations given in Table 2 helps to reduce errors in the off-diagonal 
elements. Figure 4(a) shows the experimental equilibrium spectrum and 4(b) shows the experimental 
spectrum after creating the EPR state. Since EPR state does not consist of any single quantum 
coherence, no signal is obtained (Figures 4(b)). The corresponding simulated spectrum are also 
shown on the right hand side in each figure. 

To verify the creation of the EPR state it is necessary to tomograph the complete density matrix. 
Tomography in NMR is normally carried out using spin-selective pulses obtaining a series of one- 
dimensional NMR experiments each giving a linear equation of different elements of the density 
matrix [26]. However, in the case of strongly coupled systems, all operations including tomography 
excludes the use of spin selective pulses and demands either non-selective pulses or transition selective 
pulses or both. Recently, a robust method for tomography was suggested based on two-dimensional 
Fourier spectroscopy, which utilizes only non-selective pulses [37]. This method involves: 

(i) an one- dimensional experiment for measuring diagonal elements: [G z — 10° — t], and 

(ii) a two-dimensional multiple quantum experiment for measuring all off-diagonal elements: [ti — 



The scheme involves only non-selective pulses and therefore it is not only simple and accurate but 
also applicable to strongly coupled systems. The result of measurement of diagonal elements of 
the EPR state [experiment (i)] is shown in Figure 4(c) and the corresponding simulated spectrum 
is shown in Figure 4(h). Figure 5 shows the complete pulse-sequence for creation of EPR state 
followed by measurement of the off-diagonal elements [experiment (ii)]. The resulting 2D spectrum 
of experiment is given in 4(k), which clearly shows the double quantum peaks corresponding to the 
EPR state. No zero quantum or single quantum peaks are observed. [37]. 

Since the diagonal and off-diagonal terms are measured by two different schemes, it is necessary also 
to determine the scaling between the two measurements. Normally this is achieved by comparing 
single quantum terms in the two-dimensional experiment with the spectrum obtained by direct 
detection of the density matrix [37]. However since no single quantum coherence is present in the 
EPR density matrix, we carried out two additional one-dimensional experiments: (iii)[ (tt/4)„ — t] 



(7r/2) w -G,-(7r/4)_ w -t 2 ]. 



(13) 



and (iv) [(n / A) x — t] after creation of EPR pair. Experimental spectra corresponding to (iii) and (iv) 
are shown in Figures 4(d) and 4(e) respectively, and corresponding simulated spectra are shown in 
Figures 4(i) and 4(j) respectively. The signals in 4(d) are proportional to the sum of the amplitudes 
of diagonal and double quantum coherences of the EPR state while those in 4(e) are proportional 
to the differences (when single quantum coherences are not present as in the present case). Since in 
a perfect EPR state the diagonal elements and double quantum coherences are equal, the spectrum 
of Fig 4(e) should have no signal as is evident from the simulated spectrum of Fig. 4(j). The signals 
of Fig. 4(e) compared to Fig 4(d) are measures of experimental errors, which in the present case are 
estimated to be less than 15% [31]. The complete density matrices corresponding to the theoretical 
and experimentally obtained EPR state are shown in Figures 4(1) and 4(m). 

D. Implementation of Logic gates 

Logic gates have been implemented earlier by one and two-dimensional NMR using weakly coupled 
spin-1/2 nuclei as well as quadrupolar nuclei [24,38,13,39,16]. We demonstrate here the first imple- 
mentation of a complete set of 24 one-to-one logic gates in a 2-qubit system using strongly coupled 
spin-1/2 nuclei. The system chosen for this purpose is the two phosphorus nuclei of the organometal- 
lic compound (I) shown in Fig. 6(a). The energy level diagram (Fig. 6(b)) and the equilibrium 
phosphorus spectrum of this molecule in isotropic liquid state is given in Fig. 6(c). Starting from 
equilibrium, the logic gates were implemented using sequences of transition selective pulses. The 
final populations were mapped by a small angle (10°) non-selective pulse. The spectra corresponding 
to final populations of all the 24 logic gates are given in Fig 7. The unitary transforms and pulse 
sequences for implementation of these gates are given in reference [39] with the modification that the 
r.f. power has been adjusted for given angle of flip for the two inner versus the two outer transitions. 

III. THREE-SPIN STRONGLY COUPLED SYSTEM 
A. The system and labeling of transitions 

The system chosen is the three strongly coupled protons of 3-bromo-l,2-dichlorobenzene (Figure 
8) oriented in the nematic liquid crystal ZLI-1132. The equilibrium spectrum of the system at 300 K 
obtained from DRX 500 MHz spectrometer is shown in Figure 8. There are only nine out of a total 
15 possible single quantum transitions with observable intensity, in this spin system. Construction of 
energy level diagram and labeling of transitions were performed using a Z-COSY experiment [40,41]. 



The Z-COSY spectrum along with cross-sections parallel to lu 2 axis at various transitions is given in 
Figure 9. The zero-quantum artifacts were suppressed in the Z-COSY experiment by incrementing a 
delay synchronized with ti increment [40,41]. The connectivity matrix is obtained by the MATLAB 
assisted automation 
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0-110010 
1 -1 -1 1 
-1 1 1 
1-10 0-11 
0-11001-10 
1 1 -1 
10 0-1-1100 
10 10-100 



(14) 



The constructed energy level diagram for the above connectivity matrix is shown in Figure 10. The 
ninth transition shown by dashed line belongs to the transition 011<->100, and is not connected to 
any other observed transitions. Therefore the transition did not show any connectivity to other 
transitions in Z-COSY experiment (Figure 9) and is marked as * in Figure 8. It turns out that these 
nine transitions are sufficient to carry out certain QIP operations as shown in section (B),(C) and 
(D). 



B. Preparation of Pseudopure state 

We have used the method of "POPS" to prepare a pair of pseudopure states on this three-spin 
strongly coupled system [11]. POPS requires only two population distributions: (i) Equilibrium 
populations (Fig. 11(a)) and (ii) Equilibrium populations changed by a single transition selective 
7r pulse on a given transition (Fig. 11(b)). Subtraction of (ii) from (i) yields effectively a pair of 
pseudopure states |000)(000| - |001)(001| (Fig. 11(c)). 



C. Creation of |GHZ)(GHZ| - |001)(001| state 



Entanglement between many particles is essential for most quantum communication schemes, 
including error- correct ion schemes and secret key distribution network [5]. Greenberger-Horne- 



Zeilinger (GHZ) states are three-spin entangled states of the form [44,45] 



GHZ) = -^(1000) + |111)) 



(15) 



Three particles in GHZ state exhibit one of the strangest correlations that can not be explained by 
any hidden variable theory [44,45]. A set of measurements carried out on three particles in GHZ 
state prepare the particles in a classically impossible correlated state [5]. In NMR, the GHZ state 
was first created by Laflamme et al [42]. The correlations of the GHZ state has been studied using 
NMR by Nelson et al [43]. 

The preparation of GHZ state requires preparing a pseudopure initial state, like 1 000) . However, 
since we have prepared pairs of pseudopure states as the initial state, we will be actually preparing 
a state 



This state differs from GHZ state only in the diagonal elements and therefore retains the essential 
correlations of the GHZ state. 

The GHZ state can be created from the 1 000) pseudopure state using a cascade of three transition 
selective pulses (i) pulse on the transition 8, (ii) (tt)^ pulse on the transition 4 and (iii) 

( 7r )</>3 pulse on the transition 1. The pulse sequence is, 



where pulses are to be applied from left to right. The phases 0i, 02 and 3 of the pulses can be any 
one of the 16 possible combinations, as shown in Table 3. A phase-cycle over these combinations 
helps to reduce the errors in the off-diagonal elements. 
The unitary operator for the above pulse-sequence can be written as, 



where I^J , I^J and 1^' are the single transition operators on transitions 1, 4 and 8 respectively. 
Spectrum corresponding to the diagonal part of |GHZ)(GHZ| — 1 001 ) (001 1 state is shown in Figure 
11(d). It is clearly seen that in the GHZ state (Eq. 15) transitions 1,3 and 8 have approximately 
half the intensity of the equilibrium spectrum. Transitions 5,6 and 7 appear due to POPS. However 
to confirm the creation of GHZ state a complete tomography of the created GHZ state is needed. 
The pulse sequence for preparation of POPS, creation of |GHZ)(GHZ| — |001)(001| state followed 
by tomography (using experiment (ii) of expression (13)) is given in the Figure 12(a). The 2D 



GHZ)(GHZ| - |001)(001|. 



(16) 




(17) 




(18) 



spectrum corresponding to the measurement of all off-diagonal elements of GHZ state is shown 
in Figure 12(b). Presence of only triple quantum coherence and absence of all other coherences 
confirms the creation of GHZ state. The axial peaks at zero frequency are Qi dimension is due to 
the longitudinal relaxation during the t\ period and imperfections of the 90° pulse (see caption of 
Fig. 4). 

D. Two-qubit DJ using two-dimensional NMR 

Implementation of D J- algorithm on two qubits requires three qubits including one work qubit. 
The algorithm can be described as 

\r)\s)\t)^\r)\s)\t®f(r,s)) (19) 

\r)\s) and \t) being the states of the two input qubits (hjh) and the work qubit (Jo) respectively. 
There are eight possible two-bit binary functions (f ) of which two are constant and six balanced. The 
transformations corresponding to the constant functions f\ and f 2 are respectively unity operator 
and a 7r pulse on all the transitions of the work qubit. The unitary transformations encoding the 
remaining six balanced functions f% — f$ are acheived by transition selective pulses on different 
transitions (3,4, 9(*), 6) of the work qubit as [0, 0, n, ir],[ir, rr, 0, 0],[n, 0, n, 0],[0, ir, 0, n], [ir, 0, 0, it] and 
[0,vr,7T,0]; [39]. 

The pulse sequence for two-dimensional DJ algorithm is [(|) /o ' /l,/2 — 1± — Uf — Det(t 2 )]. The transitions 
of l\ and I 2 qubits are frequency labelled during the ti period and detected during the t 2 period. 
Fourier transformation with respect to t\ and t 2 yields the desired two-dimensional spectrum [39]. 
The experimental result of the above operations on the strongly coupled 3-qubit system of 3-bromo- 
1,2-dichloro benzene is given in Figure 13. The experimental results match the expected theoretical 
results, confirming that two-dimensional DJ algorithm can also be carried out in a strongly coupled 
three spin system. 

IV. STRONGLY COUPLED 4-SPIN SYSTEM 
A. Labeling of transitions 

The system chosen is 2-chloroiodobenzene dissolved in ZLI-1132. The four aromatic protons form 
a strongly coupled 4-spin system. Equilibrium spectrum is shown in Figure 14. The transitions are 
labeled according to the descending order of their intensities. The total number of transitions for 



a 4-spin strongly coupled system is 2x4 C 4 _i = 56 of which only 30 transitions have been observed 
with sufficient intensity. There are other transitions of smaller intensity comparable to that of 13 C 
satellites. In order to avoid any interference due to 13 C satellites we have used 13 C decoupling in the 
experiment. Decoupling in t\ dimension is achieved by a single n pulse on 13 C channel in the middle 
of ti period while that in t 2 dimension is carried out by multi-pulse decoupling sequences on the 13 C 
channel. The Z-COSY spectrum obtained after 13 C decoupling for the present system is shown in 
Figure 15. The spectrum consists of more than 2000 peaks. The MATLAB analysis of the spectrum 
is carried out. Figure 16(a) shows the connectivity and labeling of the 30 transitions of Figure 14. 



B. 4-Qubit gates and pseudopure states 

Figure 16 shows implementation of gates and preparation of pairs of pseudopure states using the 
4-qubit system 2-chloroiodobenzene. The labeling scheme for the energy levels is shown in Figure 
16(a). Figure 16(b) shows the equilibrium spectrum. A C 3 -NOT gate can be implemented using 
a single n pulse on the transition 4, The spectrum corresponding to the C 3 -NOT gate obtained 
using a small angle detection pulse is shown in Figure 16(c). A pair of pseudopure states, namely, 
|1111)(1111| — 1 1 1 10) ( 1 110 1 is prepared by subtracting the spectrum 16(c) from the equilibrium 
spectrum shown in 16(b). The resultant spectrum is shown in Figure 16(d). Similarly the pair of 
pseudopure states 1 1110) (11 10 1 — 1 1010) (1010 1 (Figure 16(e)) is prepared by inverting the transition 
1 and subtracting the obtained spectrum from the equilibrium spectrum. Figure 16(f) demonstrates 
the implementation of C 2 -SWAP gate after preparing the pair of pseudopure states 1 1110) (1110| — 
|1010)(1010|. The action of C 2 -SWAP gate is to interchange the states 1 1110) and 1 1101). This is 
achieved by three transition selective pulses 



^4)^(14)^(4) 



(20) 



where the superscripts indicate the transition numbers. Since the states 1 1101) and 1 1 110) have 
almost same populations in the present case, the spectrum after applying C 2 -SWAP gate on equi- 
librium input state will not be very much different from the equilibrium spectrum. However, if one 
starts with a pair of pseudopure states 1 1110) (1110| — 1 1010) (1010 1 as the input, the output will be 
different pair of pseudopure states, 

|1110)(1110| - |1010)(1010| c2 ^ AP |1101)(1101| - 1 1010) (1010|. (21) 

The spectrum corresponding to the state 1 1101) (1 101 1 — 1 1010) (1010| obtained by applying C 2 - 
SWAP gate on 1 1110) (1110| - |1010)(1010| is shown in Figure 16(f). The pair of pseudopure state 



1 1 101) ( 1 101 1 — 1 1010) (1010| can also be prepared by inverting the transition 15 and subtracting the 
spectrum obtained from the equilibrium spectrum, as shown in Figure 16(g). The spectra in figures 
16(f) and 16(g) match fairly well, indicating good implementation of the C 2 -SWAP in Figure 16(f). 



V. CONCLUSIONS 

Increasing the number of qubits in NMR calls for the use of the dipolar couplings in oriented 
homonuclear systems which are generally strongly coupled. Spin-selective pulses are not defined in 
the case of strongly coupled systems and qubit-addressability in such a scenario is achieved through 
transition selective pulses. It has been demonstrated earlier on weakly coupled systems, that using 
only non-selective pulses and transition-selective pulses one can implement logic gates [27,39,20], and 
algorithms such as Grover's algorithm and Quantum Fourier Transform [46]. Efforts are ongoing to 
implement these algorithms in such systems as well as to realize higher qubit systems using nuclear 
spins oriented in liquid crystal matrices. 
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TABLES 



Table 1. Functions, unitary operators and r.f. pulses for one-qubit DJ algorithm. 
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Table 2. Transition numbers of expression 10 and corresponding phases of the r.f. pulses for 
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Table 3. The 16-step phase-cycle for preparing the GHZ state. The pulse-sequence is given in 

expression 17. 
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FIGURE CAPTIONS 



Figure 1. (a) Eigenstates of a strongly coupled two-spin system, (b) qubit labeling, (c) trisodiumci- 
trate yielding a strongly coupled two-spin system, and (d) the equilibrium 500MHz X H spectrum of 
(c). 

Figure 2. (a) Equilibrium spectrum of trisodiumcitrate, and spectra corresponding to various 
pseudopure states (b) 1 00) , (c) |01), (d) 1 10) and (e) |11). Numbers above the enrgy levels indicate 
populations and the binary numbers below the levels indicate the labels. Transitions are also in- 
dentified as 1, 2 , 3 and 4 in (a). The pulse sequences applied to prepare each pseudopure state is 
explained in the text. Transition selective pulses used were of length 100 ms and the gradient pulse 
was of length 1 ms and strength 10 G/cm. Each spectrum was obtained by a non-selective high 
power pulse of duration 1 fis corresponding to a flip-angle of 10°. 

Figure 3. Implementation of Deutsch-Jozsa algorithm, (a) Quantum circuit of and (b) experi- 
mental scheme for the implementation of DJ algorithm, (c) and (d) are the experimental spectra 
corresponding to the two constant functions of U\ and £/ 2 respectively, and (e) and (f) corresponding 
to the two balanced functions of U 3 and U 4 respectively. The expected spectrum for all the four 
functions are given on the right hand side. 

Figure 4. Creation and tomography of ( 1 00) + |ll))/\/2 EPR state. Spectra in (a-e) are experi- 
mental and in (f-g) are simulated corresponding to various steps of creation and tomography. (a,f) 
The equilibrium spectra, (b,g) spectra after creation EPR state (no single quantum signal observed); 
(c,h) the spectra of diagonal part measured by the pulse sequence (EPR) — G z — 10°; (d,i) spectra 
obtained by the sequence (EPR) — (ir/4) x — 10°, (e,j) spectra obtained by (EPR) — (ix/4) y — 10°. 
(k) The two dimensional spectrum to measure the off-diagonal elements. This spectrum clearly 
shows the double quantum coherence present in the EPR state. The axial peaks at zero frequency 
originate from the longitudinal relaxation of EPR state during t x period which is detected due to the 
imperfection of the 90° r.f. pulse following the t x period. The theoretical (1) and, and experimental 
(m) density matrices of EPR state. Spectra in (d,e) are used for calculating the scaling between 



diagonal and off-diagonal measurements. While plotting, the spectra shown in (d,e,i,j) are scaled up 
by a factor of 4. Pseudopure state and EPR state are created by using transition selective pulses of 
length 100 ms. An eight step cycle (shown in Table 2) was employed to minimize the errors. 

Figure 5. Pulse sequence for creation and tomography of EPR state. Numbers inside the paren- 
thesis indicate transition numbers as shown in Figure 1(d). Uepr is applied on transitions [(2), (3)] 
or [(4),(1)] as described in the expression 10 and Table 2. G\, G 2 and G 3 are field gradient pulses of 
different strengths along indirection. 

Figure 6. (a) Organometallic compound (I) in which the two Phosphorus ( 31 P) nuclei constitute a 
two spin system, (b) energy level diagram of the two spin system and (c) the equilibrium phosphorus 
spectrum recorded at 202 MHz in a magnetic field of 11.7 Tesla. 

Figure 7. Implementation of 24 one-to-one logic gates. Starting from equilibrium all the gates 
were implemented using sequences of transition selective pulses and non-selective pulses. The uni- 
tary transforms and pulse sequences for implementation of these gates are given in reference [39]. 
Gaussian shaped pulses of 100ms duration were used as selective pulses. The r.f. power has been 
calibrated for given angle of flip for the two inner versus the two outer transitions. A sine-bell shaped 
gradient was applied after implementation of each selective pulse, to kill any coherences created due 
to imperfection of pulses. The final populations were mapped by a small angle (10°) non-selective 
pulse. 

Figure 8. Equilibrium proton spectrum of 3-bromo-l,2-dichlorobenzene oriented in ZLI-1132 at 
500 MHz. The transitions are labeled from left to right. The ninth transition marked * did not show 
connectivity to other transitions (Figure 9), and belongs to the lone transition between 011 and 100 
(markedby dashed line in Figure 10). 



Figure 9. Z-COSY spectrum of oriented 3-bromo-l,2-dichlorobenzene. The equilibrium spectrum 
and the cross-sections of the Z-COSY spectrum are shown on the right-hand side. 

Figure 10. Energy level diagram of oriented 3-bromo-l,2-dichlorobenzene constructed using the 
Z-COSY spectrum shown in Figure 9. The transitions are labeled as in Figure 8. Only nine transi- 
tions are assigned. The remaining transitions are having very low intensity. 

Figure 11. Preparation of pseudopure states on the three spin strongly coupled system of Figure 
8. Energy levels, transitions and representative deviation populations (numbers inside the circles) 
are shown on the left hand side and the corresponding spectra are shown on the right hand side, 
(a) Equilibrium deviation populations and the corresponding spectrum, (b) deviation populations 
and spectrum obtained after inverting the transition 6, and (c) the deviation populations and the 
spectrum obtained by subtracting (b) from (a). The deviation populations and the spectrum in (c) 
correspond to the pair of pseudopure states: |000)(000| — |001) (001 1 . (d) Population distribution 
and spectrum corresponding to the state |GHZ)(GHZ| — |001)(001|. All spectra were recorded using 
a final small angle (10°) detection pulse to maintain linear response such that the intensities are 
proportional to the population differences of the two involved levels only. 

Figure 12. (a) Pulse sequence for creation and tomography of |GHZ)(GHZ| — |001)(001| state. 
Numbers inside the parenthesis indicate transition numbers as shown in Figure 8. G\ and G2 are 
field gradient pulses of different strengths along indirection, (b) The two-dimensional spectrum ob- 
tained by using the pulse sequence (a). Pure triple quantum coherence at iO\ + + 00$ (where, u>k 
is the frequency of the transition k) confirms the creation of GHZ state. 

Figure 13. The result of DJ algorithm on 3-bromo-l,2-dichlorobenzene desolved in ZLI-1132 for 
various functions f\ — fs- Only the expansions of the transitions of the input qubits (I\ and I2) are 
shown. Expected patterns obtained by GAMMA simulation are also shown against each spectra. 
Transitions 3,4, 9(*) and 6 are used as the four work qubit transitions. Interchange of labels 011^101 



allows us to identify transitions 4 and 9 as transitions belonging to the 3rd qubit. This does not 
affect other operations and the four 3rd qubit transitions namely 3,4,9 and 6 remain unconnected. 
However transition f now belongs to fst qubit (011<->111) along with 2 (010<->110) and transition 7 
to 2nd qubit (001^011) along with 8 (000^010). All the experiments were carried out on a Bruker 
DRX-500 spectrometer at 300K. The transition selective pulses were f .5,7.4,20 and 1.5 ms long re- 
spectively for transitions 3,4,9 and 6. The pulse power was then adjusted to make the flip angle of 
each pulse as ti. A phase cycle of (x, —x) was used to minimize the error of the 7r-pulses during 
computation. The extra peaks in / 2 (shown by " marks) in the experimental spectra were origi- 
nated due to undesired coherence transfer during computation. All the experiments were done using 
2048 t 2 and 128 ti data points. All plots are shown in magnitude mode. The resonance frequencies 
of various transitions (1,2,7 and 8) in 002 domain are schematically identified in the bottom line of 
the Figure. The tt pulses applied to various work qubit transitions are indicated for each f, with 
the transitions identified in / 2 . The same order follows for other f's. For example, / 2 = (vr, 7r, tt, ti) 
means n pulses are applied to all the transitions of work qubit and / 5 = (tt, 0, ir, 0) means n pulses 
are applied to transitions 3 and 9 and no pulses to transitions 4 and 6. 

Figure 14. One dimensional 500 MHz proton spectrum of 2-chloroiodobenzene oriented in liquid 
crystal ZLI-1132, at 300 K forming a 4-qubit system. The transitions are labeled according to de- 
scending order of their intensities. 

Figure 15. Z-COSY spectrum of the 4-spin strongly coupled system shown in Figure 13. The 
spectrum consists of more than 2000 desired peaks. 

Figure 16. (a) Labeling scheme for the states of the 4-qubit system of Figure 13. (b) Equilibrium 
spectrum obtained using a small angle (10°) pulse [the spectra of Fig. 13 was obtained using a 
90° detection pulse], (c) Spectrum corresponding to C 3 -NOT gate obtained by selective inversion of 
transition number 4 ( 1 1 1 10) <-> 1 1 1 1 1 ) ) . (d) Spectrum corresponding to the pair of pseudopure states 
|1111)(1111| — 1 1110) (1110 1 obtained by subtracting (c) from (b), named POPS(4). (e) Spectrum 
corresponding to the pair of pseudopure states |1110)(1110| — 1 1010) ( 1010 1 obtained by inverting 



the transition 1 and subtracting the obtained spectrum from the equilibrium spectrum (b), named 
POPS(l). (f) Spectrum corresponding the pair of pseudopure states |1101)(1101| - |1010)(1010| 
obtained by applying C 2 -SWAP gate on (e). The pulse sequence for C 2 -SWAP gate is given in the 
expression 20. (g) Spectrum corresponding the pair of pseudopure states 1 1101) <1 101 1 — 1 1010} < 1010 1 
obtained by inverting the transition 15 and subtracting the spectrum for the equilibrium spectrum 
(b), named POPS(15). Spectra (f) and (g) match fairly well, indicating good implementation of 
the C 2 -SWAP. All spectra were recorded using a final small angle (10°) detection pulse to maintain 
linear response such that the intensities are proportional to the population differences of the two 
involved levels only. 
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